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ABSTRACT 



We report on the multiwavelength observations of the bright, long gamma-ray burst 
GRB 110731A, by the Fermi and Swift observatories, and by the MOA and GROND 
optical telescopes. The analysis of the prompt phase reveals that GRB 110731A shares 
many features with bright Large Area Telescope bursts observed by Fermi during the 
first 3 years on-orbit: a light curve with short time variability across the whole energy 
range during the prompt phase, delayed onset of the emission above 100 MeV, extra 
power law component and temporally extended high-energy emission. In addition, this 
the first GRB for which simultaneous GeV, X-ray, and optical data are available over 
multiple epochs beginning just after the trigger time and extending for more than 800 s, 
allowing temporal and spectral analysis in different epochs that favor emission from the 
forward shock in a wind-type medium. The observed temporally extended GeV emission 
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is most likely part of the high-energy end of the afterglow emission. Both the single- 
zone pair transparency constraint for the prompt signal, and the spectral and temporal 
analysis of the forward shock afterglow emission, independently lead to an estimate of 
the bulk Lorentz factor of the jet T ~ 500 - 550. 

Subject headings: gamma rays: bursts 



1. Introduction 



Gamma-ray Bursts (GRBs) are the most powerful explosions in the Universe, initially releasing 
most of their energy in X-ray and gamma-ray on timescales lasting from a few seconds to a few 
minutes. Highly variable light curves across the energy bands in this prompt emission phase suggest 
an active central engine that drives the highly-collimated GRB jet. The prompt emission is thought 
to be emitted by internal shocks, which are produced when shells of material collide within the jet. 
The variability within this emission, as observed in the promp t light curves, is thought to be due 
to intermittent central engine activity (jRees &: Meszaroslll994l ). 



However, the details of the emission mechanism which can explain the efficiency of the in- 
ternal shocks are not und erstood. Fainter and longer-lived emiss ion following the prompt phase, 
called the GRB afterglow (jMeszaros Reeslll9971 : ISari et al.lll998l ). has been observed at lower en- 
ergies prior to Fermi, ranging from X-ray to optical and radio wavelengths. The first observations 
by t he Fermi observato ry of delayed and long-live d GeV emission relative to the prompt MeV emis- 
sion ( Abdo et al.1l2009a g) have led to speculation (Kumar Sz Barniol Duranll^Ood: 



2010 



Ghisellini et al 



2010 



De Pasquale et al. 



201 



0; 



Corsi et al.ll2010l ; iRazzaque 



Kumar fe Barniol Duran 



2010h that the af- 



terglow component may also include a significant amount of gamma-ray emission from the high 
energy tail of the synchrotron radiation of the external forward shock. A contribution from the 
fterglow to the gamma-ray flux at such early times would indicate a significa ntly earlier the onset 
of the interaction between the GRB blast wave and the circum-burst medium (Blandford McKee 



1976 ). which is thought to be the source of the afte rglow emission. The theoretical models of the 
underlying afterglow emission (IGranot fe Sarill2002l ) also constrain the physical parameters of the 
external shock such as the jet energy, bulk Lorentz factor, emission efficiency and nature of the 
surrounding medium when optical-to-GeV data are fit simultaneously in this framework. 

The Fermi observatory hosts two instruments, the Lar ge Area Telescope (L AT) which covers 
an energy range from 20 MeV to up to more than 300 GeV (jAtwood et al.ll2009l ) an d the Gamma 



ray B urst Monitor (GBM) which is sensitive at lower energies, from 8 keV to 40 MeV (IMeegan et al 
2009). Together, the LAT and GBM are capable of measuring the spectral parameters of GRBs 



across seven decades in energy. During the first 3 years on-orbit, the Fermi LAT has detected more 
than 30 GRBs with high significance, but only four of these have benefited from simultaneous detec- 
tions with the Swift Burst Alert Telescope (BAT), allowing for prompt follow-up wit h the narrow 



field X-Ray Telescope (XRT) and Ultraviolet Optical Telescope (UVOT) on Swift (jGehrels et al 
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2004). The first such burst was GRB 090510, a very bright short GRB for which the afterglow was 
observed contem poraneously by Fermi LA T and by Swift XRT and UVOT starting just 100 s after 
the burst trigger (|De Pasquale et alj|2010l ). Multiwavelength data from GRB 090510 seem to favor 
an a fterglow interpretation as their origin (|De Pasquale et al.l |2010| ; ICorsi et al.l 120101 ; iRazzaque 
2010) . The second such burst was GRB 100728A, a very long burst for which Fermi revealed sig- 
nificant temporally extended emission out to 850 s post trigger, and Swift detected a series of strong 
X-ray flares in the XRT light curves. An internal shock scenario seems to reproduce well both the 
prompt emission an d the later X-ray flares as well as temporally ext ended high-energ y emission 
from GRB 100728A dAbdo et alJboilh . The third was GRB 110625A (|Tam et allboij ) which is a 
burst that triggered both BAT and GBM but which was outside of the LAT field-of-view during 
the prompt phase. However, thanks to the Fermi Autonomous Repoint Request (ARR) there was 
a joint XRT and LAT detection of the temporally extended emission. 

In this paper we report on the analysis of the bright and long GRB 110731 A, the fourth 
burst to benefit from joint Fermi and Swift observations, and one with the most comprehensive 
multiwavelength data of any LAT-detected GRB to date. Fortuitously, the burst was within the 
LAT field of view at the trigger time and caused an ARR. BAT also triggered on the burst and 
Swift was immediately repointed for XRT and UVOT observations. In addition, the burst was also 
observed by ground-based observatories, including the Microlensing Observations in Astrophysics 
(MOA) telescope at early and intermediate times, and the Gamma-ray Burst Optical/Near- Infrared 
Detector (GROND) at late times. Multiwavelength observations hence cover both the prompt phase 
and the temporally extended emission of GRB 110731A, from optical to GeV energies. 

The paper is organized as follows: Section [2] includes a comprehensive description of the 
observations made by the various instruments that detected the burst's prompt and temporally 
extended emission, Section [3] gives the details of the data reduction and analysis for both the 
prompt emission and temporally extended emission for all telescopes, Section U] reports the results 
of the prompt and temporally extended emission data analysis in the multiwavelength context and 
Section [5] provides the discussion and interpretation of both the prompt and temporally extended 
emission results, and finally we draw our conclusions in Section [6l 

Throughout the paper, times t = T — Tq are given relative to the GBM time of trigger To, and 
the afterglow convention for the energy flux F v t oc v~@t~ a has been followed, where the energy 
index f3 is related to the differential photon index r = /3 + 1. The phenomenology of the burst is 
presented in the reference frame of the observer, unless otherwise stated. All the quoted errors are 
given at the 68% confidence level for one parameter of interest. 



2. Observations 



On 2011 July 31 at 11:09:29.94 (UT), GBM triggered on GRB 110731A which, due to the high 
peak flux of this burst, caused an ARR. GRB 110731A was already well within the LAT field of 
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view, being only ~ 3.3° off axis so that the repointing had little impact on the prompt emission 
phase observations. The maneuver placed the spacecraft in pointing mode for 2.5 hours after the 
burst, allowing continuous LAT observation of the burst from the initial time of trigger until the 
first Fermi passage into the South Atlantic Anomaly (SAA) at 1400 s. The ARR continued for 
another 90 minutes after Fermi had exited the SAA at 3150 s, with the burst being well within 
the LAT field of view from 4000 s to 7400 s. No significant signal from GRB 110731 A was found 
at these late times. 



The best LAT localization of GRB 1 10731A is RA=18 1 | 1 41 m 00 s , DEO 
with a 68% confidence error radius of 0.2° (jAbdo et al.ll2009bl ) 



-28°31'00" (J2000), 



GRB110731A triggered the Swift BAT at 11:09:30.45 UT. Swift slewed immediately to the 
burst, and its narrow-field instruments, XRT and UVOT, began observations 56 s after the BAT 
trigger. An accurate afterglo w position w as rapidly determined by the UVOT as RA=18 h 42 m 00.99 s , 
DEC=-28°32'13.8" ( J2000. lOates et alJEfflllh . with an error radius of 0.5 arcsec (90% confidence). 



XRT observations started while the spacecraft was settling at the end of the initial slew. The 
XRT began collecting data in Window Timing (WT) mode, as the source was bright (~100 cts s _1 ), 
and automatically switched to Photon Counting (PC) mode when the count rate from the source 
decreased to < 2 ctss -1 . Follow-up X-ray observations occurred during the following 24 days for 
a total net exposure of 600 s in WT mode and 75 ks in PC mode. 

The UVOT took a short exposure with the v filter during the settling phase. This exposure 
was followed by a 'finding chart' exposure with the White filter lasting 147 s. UVOT then began 
its usual procedur e of cycling through its 3 visible filter s (v, b, and u) and 3 UV filters (uvwl, 
uvm2, and uvw2) (jPoole et al.l 120081 : iBreeveld et al.l 12010 ) . The optical afterglow was detected in 
the White, u, b, and v filters, but not in the UV filters. The lack of detection in the UV filters is 
consistent with the measured redshift of z = 2.83 (jTanvir et al.ll201ll ). 



MOA observatio ns began 3.3 min after the Swift trigger for GRB 110731A (GCN 12242, 



Tristram et al. 



2011 



). Using a 61 cm Boiler & Chivens telescope at the Mt. John University 
Observatory in New Zealand, I and V band images with 60 s exposures followed by 120 s exposures 
until 12:56 UT (105 min after the trigger) were obtained. The total numbers of I and V images are 
39 and 35, respectively; however due to the difficulty of photometry in the crowded field, we used 
only the 30 I and 19 V data points reported in Tab. [73 



After a weather-induced delay, the seven-color imager GROND (jGreiner et al.ll2008l ) mounted 
on the 2.2 m MPG/ESO telescope at La Silla Observatory, Chile, observed GRB 110731A at a 
mean time of 2.74 days after the trigger. Two 30-minutes observation blocks were obtained which 
yielded an integration time of 4500 s in g'r'i'z' and 3600 s in JHK. The mean seeing during the 
observations was l'/3. 
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3. Data Analysis 
3.1. Fermi 

For the time-integrated and time-resolved spectral analysis of the prompt phase, both GBM 
and LAT data were used. The Fermi LAT and GBM data may be retrieved from the Fermi Science 
Support Center archives 0, 1^1 . 



The GBM detectors were selected in the same fashion as outlined in (jAbdo et al.1 l2009bl : 



Gruber et al.ll201ll ; iGoldstein et al.ll2012l ): we used the sodium iodide (Nal) detector s and 3, and 



bism uth germanate (BGO) detector 0. We also used Time Tagged Events TTE data (jMeegan et al 



2009) for our spectral analysis with a temporal resolution of 64 ms, in the 8 keV to 40 MeV energy 
range, excluding the range around the Nal K-edge at 33.17 keV. 

For the LAT, we first extracted 'P7TRANSIENT'-class data from a circular region centered on 
the burst position with energy-depend ent radius equal to a 95% containment of the point-spread 
function (PSF), see lAbdo et al.l (|2009bl ) for details. To greatly reduce the numbers of gamma rays 
from the Earth limb w e selected ev e nts wit h zenith angles less than 100°. We then followed the 
procedure described in lAbdo et al.1 (|2009bl ) to estimate the residu al background and considered 
front- and back-converting events separately (| At wood et al.l |2009| ) . The results presented here 
were obtained using the Fermi ScienceTools-v9r25pl and P7TRANSIENT_V6 instrument response 
functions (IRFs). 

Besides the 'P7TRANSIENT' data, we also extracted events u sing the so-called LAT low 
energy (LLE, E> 10 MeV) events selection criteria (jPelassa et al.ll2010i ): This selection keeps events 
that pass the GAMMA filter, have a reconstructed track in the tracker pointing to a sky location that 
is roughly compatible with the GRB position (jAtwood et al.ll2009i ). By retaining very low energy 
events, this selection provides high statistics light curves that are useful for temporal analysis, see 
Section 14.11 We did not use LLE data for the spectral analysis since LLE and 'P7TRANSIENT' 
data gave consistent results. 

The LA T /GBM joint spectral fits were performed with the software package XSPEC ver- 
sion 12.7.0e ( Arnaudlll996l ). and cross-checked with rmfit version 4.0rcl [3 (IKaneko et al.l 12009 : 
Abdo et al.ll2009bh . GB M Response Matrice s vl.8 were used with both tools. For the fitting pro- 
cedure, the PG-statistic f Arnau d et al.ll20 11) (S in the following) was used within XSPEC, while 
the Castor statistic (jAckermann et al,ll201ll ) was used for rmfit. 



For each time interval of interest, we compared the fit of several models that are listed below 



77 http: / /heasarc. gsfc.nasa.gov/FTP/fermi / data/gbm / triggers / 
78 http://fermi. gsfc.nasa.gov/ssc/data/access/ 



79 



rmfit for GBM and LAT analysis was developed by the GBM Team and is publicly available at 
fermi.gsfc.nasa.gov/ssc/data/analysis/ 
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(differential fluxes are in ph/cm 2 /keV/s): 



(i) 



a power law with an exponential cutoff (Comptonized model, hereafter COMP), whose dif- 
ferential photon flux is described by the equation: 



COMP(£) = A 



E 



E, 



ref 



(i) 



where A is the normalization amplitude, T a is the photon index, E re f is the reference energy 
fixed at 1 MeV and Eq is the cutoff energy. 



(ii) 



a Band function, hereafter BAND, as defined in ([Band et al.lll993l ). where A is the normal- 
ization amplitude, T a and are the low and high energy power law indices and Eq is the 
cutoff energy: 



BAND(£) 



A dookev) ex P(~i|' 



.4 



(r Q -!>)£() 

100 keV 



, 100 keV 



E < (T a - Tfs)Eo 

) r/3 ex P (r^-r Q ) E>(r a -T f) )E 



(2) 



(iii) 



a COMP model plus a power law (COMP+PL): 

COMP+PL(£) = COMP(£) + B x 



E 



(iv) 



100 keV 

where B is the normalization and T 7 the power law photon index, 
a BAND model plus a power law (BAND+PL): 

BAND+PL (£) = BAND(£) + B x ' 



(3) 



100 keV 



(4) 



(v) a model with a sum of two independent Comptonized components (COMP+COMP) that in 
addition to the extra PL component provides also a second high-energy cutoff. 

We remind the reader that the commonly referenced peak energy -Epeak of the vF v spectrum can 
be computed as: 

Epeak = (2 + T a )E in the COMP model. (5) 

The COMP (i) and BAND (ii) model s reported above were historically foun d to give very good 
empirical descriptions of GRB spectra (jBand et al.l Il993l ; iKaneko et al.l 120061 ) . with the COMP 
model having one less free parameter and providing a sharper decrease at high energies (usually in 
the ~MeV range) than the BAND model. In order to check whether a different spectral component 
(meaning a different physical process and potentially a different location of the emission) could be 
responsible for the GeV emission observed in the LAT data, we added a power law component to 
models (i) and (ii), resulting the new models COMP+PL (iii) and BAND+PL (iv). Additionally, 
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as some theoretical models predict that the high-energy GeV emission must be cut off owing to 
gamma-gamma pair production opacity considerations, see Section 15.14 we tested this hypothesis 
by considering the COMP+COMP model (v), that matches these characteristics while keeping a 
reasonable number of free parameters. 

For each time interval we determined the best-fit model that provides a good description 
of the data with a minim al set of parameters, following the method described in Section 4.1 



of lAckermann et al.l (|201ll ). In particular, we used the likelihood ratio test (LRT) to derive the 
significance of the improvement of the fit when comparing a simpler model (the null hypothesis) 
with a more complex model (the alternative hypothesis). Using the PG-stat value, S, of the statis- 
tic and defining AS as the difference between the values of S obtained with the two models, the 
LRT gives the probability P(AS) that the observed AS has been obtained because of statistical 
fluctuations, on the assumption that the null hypothesis is the true model. Thus, if P(AS) is low 
the alternative hypothesis is to be preferred. We sampled the distribution of AS* via 10 millions 
Monte Carlo realizations of the burst spectrum with XSPEC. With such large statistics, we were 
able to compare pairs of models for probabilities down to P(AS) ~ 1 x 10~ 7 . Such simulations 
cannot account for systematic effects, for example due to the uncertainties in the responses of the 
instruments. Although a 10 -4 probability would be formally very significant, we adopted a conser- 
vative threshold P t h = 1 X 10 -5 , and we preferred the alternative models over the null hypothesis 
if P(AS) < P th . 

Since the count flu ence of this burst is very high, we used the effective area correction fac- 



tors (lAbdo et al.ll2009al ) to account for possible calibration issues between the different detectors. 
We used the time-integrated spectrum (interval PI, see Section I4.1.2P and the BAND model to 
determine these factors, which were 0.85 for the BGO detectors and 0.96 for the Nal detectors, 
using the LAT spectrum as a reference. We applied these effective area corrections throughout the 
spectral analysis. 

We explore the temporally extended emission from the burst in the LAT data running an 
unbinned likelihood analysis with the Fermi ScienceTools-v9r25pl and P7TRANSIEN T_V6 IRFs. 



Again, we estimated the residual particle background following the procedure described in lAbdo et al 



(|2009bl ) for each sub-data selection. The model prepared for the likelihood analysis was then com- 
posed of a single power law component for the burst itself, the particle background template and 
the template model 'gal_2yearp7v6_v0. fits' for the Galactic diffuse emission for which the nor- 
malization was held fixed during the fit. 

In order to test the significance of the detection/light curve of the GRB at a particular time, 
without being dependent upon a subjective choice of time interval, we determined interval bound- 
aries through the use of an algorithm that extends the time intervals until the likelihood test statistic 



80 Fermi background models are available from the FSSC web site at 
http://fermi.gsfc.nasa.gov/ssc/data/access/lat/BackgroundModels.html 
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is 18 in each bin. If a time interval exceeds the good time interval, an upper limit is reported and 
the computation steps to the next good time interval. The minimum number of events required in 
each interval is 8, to guarantee a reasonable number of degrees of freedom in the fit. Note that this 
method also optimizes the detection probability and and provides the m aximum number of time 
intervals on which the analysis can be run (more details may be found in lLott et al.l (|2012l )). 



In order to fit the spectral energy densities (SEDs) using data from the different instruments 
(see Section [4.2.21 below) . it was necessary to combine and adjust the LAT time intervals to have 
sufficient counts per bin and to match the observed features of the XRT light curve (see Section 14.21 
below). We defined the new time intervals as follows: [8.3, 11.5] s (I, hereafter), [11.5, 55.0] s (II), 
[55.0, 227.0] s (III) and [227.0, 853.9] s (IV). We fit the LAT again over these time intervals using 
a slightly different background model in which the particle background and th e Galactic diffuse 
emission are both estimated using the procedure described in lAbdo et al. ( 2009b ) so that we could 
extract separately the signal and background data to be used as input to XSPEC for the joint fits 
of the SEDs. 



3.2. Swift 



We retrieved the Swift data from the HEASARC archival and processed them with the 
standard Swift analysis software (v3.8) included in the NASA's HEASARC software (HEASOFT, 
ver. 6.11) and the relevant calibration files. 

We extracted BAT mask-weighted l ight curves and spectr a in the nominal 15-150 keV energy 
range following the standard procedure (jSakamoto et al.ll2008l ). The BAT data were not used for 
the spectral analysis of the prompt phase because the GBM data alone constrain r a very well and 
because the cross-calibration between BAT and GBM is still not well understood and the subject 
is beyond the scope of this paper. 

We extracted the XRT light curves and spectra in the nominal 0.3-10 keV energy range by 
applying standard screening criteria. All the XRT data products presented here are background 



subtracted and corrected for PSF losses, vignetting effects and exposure variations (see lEvans et al 



2007 



2009J). 



UVOT photometric measurements were complicated by the crowded field. We obtained the 



source count rates fr om a circular source extraction region with a radius of 5" (jPoole et al 



2008 



Breeveld et al.1 1201G ) . We estimated the background from nearby circular regions with radii of 
20", whithin which field sources were masked out. We also used Swift late-time observations to 
estimate the residual contribution of nearby objects, and to refine the afterglow photometry. For 
this reason, the last data point used in the light curve is at 1000 s. In order to better constrain 



http://heasarc.gsfc.nasa.gov/docs/swift/archive/ 
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the optical temporal decay we created a single li ght curve (see Fig . H]) from all the UVOT filters 
by re normalizing each ligh t curve to the w-band (jOates et al.ll2009l ). using flux conversion factors 



from iBreeveld et al.1 (|201ll ). 



3.3. MOA 

We performed the MOA data analysis using aperture photometry via the SExtractor pack- 



age (jBertin &: Arnoutsl ll996) to estimate the instrumental magnitudes of the objects in each im- 
age. We then compared these values to the late time GROND data of the same field for zero point 
determination of each of the MOA images. Since the MOA data were obtained with Bessell's / and 
V band filters, which have transmission curves similar to those of Johnson-Cousin's / and Fband 
filters in combination with the CCD q uantum effici e ncy c urve, and GROND uses the g'z'r'i' filter 



system, we used the conversion table in ISmith et al.1 J2OO2I) to obtain a GROND equivalent / and V 



band magnitudes with which to perform this comparison. In order to reduce the systematic error 
in the zero point determination, we selected only the top 30th percentile of brightest stars in the 
MOA images. Finally, we produced light curves for the same bright stellar objects in each of the 
MOA images for a final relative calibration. The resulting median Am variations in these stellar 
light curves were used as additional corrections to the afterglow light curve that account for any 
errors in the zero point determination of the individual images. We estimated the systematic error 
in the zero point calibration of each image by measuring the standard deviation of the difference 
between the MOA and GROND equivalent / and V band magnitudes in each image. We then 
summed the statistical errors of the measured fluxes as returned by SExtractor in quadrature with 
this systematic error to obtain the error estimate in the final MOA flux density data. MOA results 
are given in Tab. and as fluxes reported on the multiwavelength light curve shown in Fig. 01 



3.4. GROND 



We an alyzed the GR OND data analysis with a custom routine as described in iKriihler et al 



(|2008l ) and lYoldas et al.1 (120081 ). using SExtractor for background subtraction, and masking out 
bright sources. At the position of the afterglow, a faint source was visible in r'i'z' , but photometry 
was hindered by multiple nearby stars in the crowded field. Therefore, we obtained obtained mea- 
surements for a second epoch on 2011 September 25, with identical exposure time, under improved 
conditions, to create a template image for image subtraction. We performed the image subtrac- 
tion using HO TP A NT$K We used 60 different parameter settings in determining the Gaussian 
PSF kernel for the subtraction routine for each band, and chose the combination of parameters 
that resulted in the best subtraction of nearby stars near the afterglow position, as measured by 



http://www.astro.washington.edu/users/becker/hotpants.html 
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the noise in the residual image at the afterglow position. In the residual image, the afterglow is 
strongly detected in r', still well-detected in i'z', only faintly detected in g'JH and undetected 
in K. We calibrated the magnitudes of stars in the field against an SDSS standard star field at 
similar RA observed just before the first epoch observations under photometric conditions. We per- 
formed the photometry using seeing-matched aperture photometry on the subtracted images with 
MIDAlQ We estimated the errors on the fluxes as the sum in quadrature of the calibration error, 
the statistical error of the detection, and the noise error of the image subtraction as determined by 
HOTPANTS. GROND results are given in Tab. El and as fluxes reported on the multiwavelength 
light curve graph shown in Fig. |U 



Results 



4.1. Prompt phase 



4-1.1. Light curves and Timing results 

In Fig. [H we show the GBM and LAT light curves of the GRB prompt emission phase in 
several energy bands, from 8 keV to above 1 GeV. The light curves show a complex multi-peaked 
structure and have two interesting features: (i) the LAT emission at > 10 MeV is slightly delayed 
(~ 2.5 s) with respect to the GBM light curves, (ii) a peak with high count rate is also present at 
5.5 s in the LAT data that is also present in the Nal and BGO light curves. 



Detailed analysis of the GBM data results in a T90 gbm duration (|Kouveliotou et al.lll993l ) in 
the 50 keV to 300 keV energy range of 7.3±0.3 s, with a start time defined by T^gbm = 0.25±0.1 s 
and an end time at Tgs^BM = 7.6 ± 0.3 s. A similar detection and duration analysis of the LLE 
light curve (> 10 MeV) demonstrated that the LAT prompt phase detection starts later than 



14.3~f 7 D s. These 



observed in the GBM, at Tos,lle = 2.5_ ' 6 s, and lasts longer with Tgo,LAT 
results for the LLE event selection are consistent with the findings for 'P7TRANSIENT'-class 
events (> 1 00 MeV) using an es timation of the total background following the procedure described 
in detail in Ubdo et all (j2009bl ). Indeed, for the 'P7TRANSIENT' events we found a comparable 



Tt)5,LAT 
s. 



3+0.3 
5 -0.2 



s. This analysis also revealed temporally extended emission up to Tgs^AT 



1901™ 



The LLE light curve was investigated using a Bayesian blocks algorithm ( Scarglefl998 : Scargle et al 



2012J) to determine intervals over which the photon arrival rate has no statistically significant vari- 
ations. By requiring a large statistical significance for the rate variations (to be less sensitive to 
background noise) four time intervals were found at [0, 2.44, 5.44, 6.06, 8.52] s and are highlighted 
in Fig. [TJ The upper boundary of the first time interval is at 2.44 s which is consistent with the 
previous estimate of the first signal significant detection in the LAT as reported above with other 



http: / /www.eso.org/ sci/software / esomidas / 
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Fig. 1. — GBM and LAT light curves for the gamma-ray emission of GRB 110731A. The data 
from the GBM Nal detectors were divided into soft (8-14.3 keV) and hard (14.3-260 keV) bands to 
reveal similarities between the light curves at the lowest energies and for the LAT data. The first 
four light curves are back ground-subtracted and have 0.1 s time binning. The fourth panel shows 
the LAT LLE light curve dPelassa et aliboioh . The fifth panels shows LAT T7TRANSIENT'-class 
events light curves for energies > 100 MeV, with 0.5 s time binning. The sixth panel shows the 
energy and arrival time of LAT 'P7TRANSIENT'-class events above 100 MeV The vertical lines 
indicate the boundaries of the intervals a, b, c, d, used for the time-resolved spectral analysis for 
boundaries [0.00, 2.44, 5.44, 6.06, 8.52] s. The insets show the counts for each data set, binned 
using these intervals, to illustrate the numbers of counts considered in each spectral fit. 
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techniques. The four time intervals will be used in the next Section for the time-resolved spectral 
analysis. 

The flux peak at 5.5 s is observed across the whole energ y spectrum so we exam ined it by 



calculating the cross-correlation functions (CCFs), as defined in lEdelson fc Krolikl (|1988l ). between 
energy bands to quantify the simultaneity of the emission. We calculated the CCF between the 
low energy Nal and the LLE light curves, first considering the full prompt phase from — 1 s to 
8.52 s: the CCF has a local maximum at 0.3±gf s, consistent with no lag. However, a second CCF 
maximum occurs at 5.2 s, possibly due to overlapping pulses and making the lag measurement 
subject to systematic uncertainties. The CCF thus suggests that the peak is indeed statistically 
coincident in time across the whole energy range, although it is less prominent at the very lowest 
energies in the Nal detectors than at high energies (the LAT events above 100 MeV). 

We further characterized the variability of the emission by deriving a typical variability time 
sc ale for the burst emission following the light curve pulse deconvolution technique described 
in lBhat et al.l (|2012l ). We used the summed light curves from the 4 brightest Nal detectors, both 
the BGO detectors, as well as the LLE data. The light curves were sub-divided into various en- 
ergy bands, when possible, to estimate the energy dependence of the full width at half maximum 
(FWHM) of the pulses. The median value of the FWHM of the pulses in the GBM Nal and BGO 
data over the Tgo,GBM interval shows a weak dependence on energy, decreasing from 0.35 ± 0.02 s 
at ~ 18 keV to 0.24 ± 0.01 s at 2 MeV with an average value 0.28 ± 0.02 s. The median value of 
the FWHM in the LLE data during the T 90 ,gbm interval is larger, 0.92 ± 0.15 s at ~ 17 MeV. In a 
shorter time interval, [3.0 s, 7.6 s], that excludes the Tos^at time, the median value for the FWHM 
is 0.43 ± 0.03 s with a minimum of 0.147 ± 0.003 s, in the 8 keV-1 MeV range. 

For each LAT 'P7TRANSIENT'-class event, we estimated the probability of its being associ- 
ated with the GRB using the gtsrcprob Fermi Science Tool. The probability computation takes into 
account the spectral and spatial distributions of all of the components in the source model, con- 
volved with the response of the LAT as well as the exposure (all convolved with the e ffective PSF). 



The values of the model parameters are found via a maximum likelihood analysis (jMattox et al 



19961 ). and the probability of a particular event being attributed to a particular source component is 
proportional to the predicted counts density for the event by that component. The highest photon 
energy during the prompt phase is a 2.0 GeV event at 8.27 s, having a probability of 10~ 6 to be 
associated with the background. During the temporally extended emission phase, we note that a 
3.4 GeV event at 435.96 s has a probability of ~ 10~ 3 to be associated with the background. 



4-1-2. Spectral analysis 
We analyzed the burst emission spectrum over a number of time intervals: 



- the 4 intervals a, b, c, d which were determined by the Bayesian blocks analysis in Section f4. 1.11 
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- interval PI (a + b + c + d or [0, 8.52 s]) corresponds to the entire prompt emission phase. 

- interval P2 ([3.0 s, 7.6 s]) spans the time range To5,lat to I^gbm in which we observe the 
maximum flux in both instruments. 

We first considered the time interval PI for the time-integrated spectral studies. We found 
that BAND is the preferred model. Both COMP+PL and BAND+PL provided only limited im- 
provements in the fit. Comparing each one separately with BAND, we obtained a null hypothesis 
probability P(AS) ~ 1 x 1CT 4 for both. The best-fit parameters for BAND were T a = -0.89 j]^, 
r^g = -2.32jl^o|, E = 324+25 keV with a corresponding fluence in the 10 keV - 10 GeV energy 
band of F = (4.56 db 0.05) x 10 -5 erg cm -2 , calculated in the rest frame. 

To further investigate the significance of the additional power law component, we performed a 
time-integrated spectral analysis on interval P2. Similarly to interval PI, the BAND model fitted 
the data reasonably well. This time, however, BAND+PL gave a large improvement and smaller 
residuals with respect to BAND, with AS = 35.2 and a corresponding null hypothesis probability 
of P(AS) < 1 x 10~ 7 . Thus, the power law component is required to account properly for the high- 
energy part of the spectrum. Moreover, the COMP+COMP model provided an even better fit for 
the same number of degrees of freedom, with AS = 15.5 when compared to BAND+PL. Using 
the latter as null hypothesis we obtain P(AS) ~ 3 x 10~ 5 . This latter result strongly suggests 
the presence of a cutoff in the high-energy part of the spectrum, although the null hypothesis 
probability is not formally below P th . A COMP+BAND or BAND+BAND model did not improve 
the fit over the COMP+COMP model. Furthermore the power law slope of the high-energy BAND 
component could not be constrained due to limited statistics. The best-fit parameters of this time 
interval are given in Tab. [TJ and the count spectrum corresponding to BAND+PL is shown in 
Fig. El 

We performed the time-resolved analysis in the 4 intervals a, b, c, d determined by the Bayesian 
blocks analysis (see Section [4. ITT]) : a Bayesian blocks algorithm ensures that the flux has no statis- 
tically significant variation over each interval, and hence provides a useful binning for studying the 
spectral evolution of the burst, in particular in the MeV-to-GeV range because the algorithm was 
run on the LLE light curve. We investigated in particular the presence of the power law component 
and of a high-energy cutoff. The best-fit spectral parameters for the preferred models are reported 
in Tab. El and the best models both for the time-resolved and time-integrated analyses are shown 
in Fig. El The power law component is not detected above the probability threshold in any of the 
intervals, likely due to the limited statistics. The results suggest, though, that a below-threshold 
power-law component may be present starting after the second interval. There is also marginal 
evidence of a high-energy cutoff in the second interval. Here we summarize the results from the 
time-resolved spectral studies: 

- interval a ([0, 2.44 s]): the only model that was constrained by the data is COMP. 
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Fig. 2.— Joint spectral fitting of GBM and LAT data for the time interval P2, [3.0 s, 7.6 s]. The 
top panel shows the count spectra (points) and best-fit BAND+PL model (lines). The lower panel 
shows the residuals. 
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Table 1: Best-fit parameters for all the models for the spectrum obtained in the interval P2, [3.0 s, 
7.6 s]. The reference energy is fixed to 1 MeV. The fluence is evaluated for the 10 keV-10 GeV 
range. 

- interval b ([2.44 s, 5.44 s]): Adding a PL to BAND gave an improvement AS = 19.5. Adopting 
BAND as null hypothesis, the corresponding probability is P(AS) ~ 2 x 10 -5 , very close 
to the significance threshold. COMP+COMP marginally improved the fit with respect to 
BAND+PL, giving AS = 8.1. Adopting the latter model as null hypothesis, we obtain 
P(AS) ~ 3 x 10" 4 . 

- interval c ([5.44 s, 6.06 s]): BAND, BAND+PL and COMP+COMP gave very similar values 
for S. Thus, from a statistical point of view, the more complex models are not favored over 
BAND. However, it is noteworthy that the best-fit parameters of the complex models gave 
results which are consistent with what we found for interval b. 

- interval d ([6.06 s, 8.52 s]) : BAND+PL fit the data slightly better than BAND alone, with 
AS = 13.7. Adopting BAND as null hypothesis, we obtain P(AS) ~ 3 x 10~ 4 . The high- 
energy cutoff of the COMP+COMP model is not well constrained. 



4.2. Afterglow modeling 

4-2.1. Multiwavelength afterglow light curves 

The multiwavelength light curves for the 7 instruments (LAT, GBM, XRT, BAT, UVOT, MOA 
and GROND) that observed GRB 110731A are shown in Fig. [U and the corresponding data points 



-19- 




10 1 2 1 3 1 4 1 5 1 6 1 7 

Energy (keV) 



Fig. 3. — Top: The best-fit BAND+PL model for the time-integrated interval P2 plotted as a 
vF v spectrum. The two components are plotted separately as the dotted lines, and the sum of 
the components, representing the overall spectrum, is plotted as the heavy line. The ±1 a ranges 
derived from the errors on the fit parameters are also shown. Bottom: The uF v model spectra (and 
±1 a error contours) plotted for each of the time bins considered in the time-resolved spectroscopy. 
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Time interval from T [s] a (0-2.44) b (2.44-5.44) c (5.44-6.06) d (6.06-8.52) 
Best model COMP BAND BAND BAND 

E [keV] I88tjj 285±|g 683±f|g 446±?| 

r n no+0.05 n c/i+0.05 i 1 c+0.05 n o fi +0.06 

1 Q -U.y2_ 005 -U.b4_ 05 -1.15_ 006 -U.Sb_ 006 

r o O/I+0.04 9 1s +0.05 901 +0.04 

L P - -^-34„ 04 -2.18_ 06 -2.31_ 0Q5 

Fluence [10- 5 ergcm- 2 ] 0.58±g;g| 2.05±g;gJ 0.59i°;°! 1.10±g;g| 

PG-stat (DOF) 

BAND (354) - 417.4 365.3 389.1 

COMP (353) 378.6 

BAND+PL (352) - 397.9 363.2 375.4 

COMP+PL (353) - 399.7 365.3 380.2 

COMP+COMP (352) - 389.8 360.2 377.7 

Table 2: Summary of GBM/LAT joint spectral fitting by best model in 4 time intervals. The 
fluences are evaluated for the range covered by both instruments is 10 keV - 10 GeV. 

are reported in Tab. H] through Tab. [8] in appendix. 

The LAT light curve in Fig. [4] shows that the peak of the flux density is in the interval 
[5.47 s,5.67 s], and that after 5.67 s, the flux decays smoothly during the whole temporally extended 
emission, following a power law with a fitted index a = 1.55 ± 0.20. The last time interval with a 
clear detection in the LAT data spans from 227.0 s to 853.9 s. 

The BAT light curve shows many consequent adjacent peaks, from the trigger time to ~8.5 s. 
Faint, spectrally soft emission is detected out to ~30 s, and its temporal decay can be described 
by a power law with index obat=2.3±0.3. By shifting the reference time to the peak of the 
observed LAT emission at t p k=5.5 s, i.e. shifting the reference to the onset of the forward shock, 
the decay slope changes to aBAT,t pfe =l-25±0.15. A rebrightening is visible between 50 s and 80 s, 
in coincidence with the first X-ray flare as seen on the multiwavelength light curve (Fig. U|) . 

The X-ray light curve shows two X-ray flares, peaking at 70 s and 110 s respectively. By 
excluding the interval of significant variability (from 56 s to 150 s) from the temporal fit, the light 
curve can be described by a simple power law with decay slope ax = 1.189 ±0.007 (x 2 =250 for 228 
d.o.f.). This model, though acceptable at a statistical level, systematically overestimates the late 
time flux densities. A broken power-law yields a significantly better result (x 2 = 217 for 226 d.o.f.) 
and shows no systematic trend in the residuals. The best-fit parameters are axi = 1.10 ± 0.02, 
a x ,2 = 1-32 ± 0.03 and t bk = 4.6±^J ks. 
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Fig. 4. — Multiwavelength light curves of GRB 110731A observations: vertical dotted lines define 
the time boundaries of the SEDs (indices I, II, III, IV and V) studied in Section 14.2.21 The power- 
law index of the LAT light curve is «lat = 1-55 ± 0.20. The XRT light curve follows a broken 
power law with ax,i=l-10±0.02, ax,2=l-32±0.03 and a break at t bk =4.6tH ks - The UVOT light 
curve is well fit by a single power law with decay slope a op t=l-37±0.03. 



-22- 



For the optical emission, the light curve we obtained using UVOT and MOA data altogether is 
very well fit by a single power law with decay slope a op t = 1.37 ± 0.03, as shown in Fig. 21 For the 
UVOT data, we created a single filter light curve, the White light curve, by co-adding normalized 
version of the light curves for the 7 UVOT filters. The multiwavelength light curve in Fig. 0] shows 
just UVOT White and v data for clarity, as well as MOA data for the / and V bands. Data points 
for individual and co-added exposures are reported in Tab. [5] and Tab. [6l 



4.2.2. Broadband SED 



We evaluated the SED at 5 epochs: 9 s (I), 16 s (II), 100 s (III), 550 s (IV), and 2.74 d (V). 
The SED times were selected in order to maximize the data coverage and minimize interpolation. 
For each epoch and instrument, a spectrum was extracted using the corresponding time interval, 
and scaled to the actual count rate at each time of interest. The time intervals were chosen in order 
to achieve a significant detection in the LAT energy band (epochs I-III) , sufficient statistics in the 
X-ray spectrum (epoch IV) and to just match GROND observations for epoch V. Optical fluxes 
were interpolated, when necessary, by using the best-fit model for the light curve. At the GRB 
redshift of z = 2.83 the Lyman series absorption is redshifted to the observer frame wavelength 
range 3500-4700 A. This mainly affects the UVOT b, u, and the GROND g 1 data points, which 
were therefore excluded from the spectral fits. 

Each SED was fit with an absorbed power law or an absorbed broken power law, where the 
spectr al slopes 0-\ and 0i were ti ed to 0i = 0\— 0.5, as predicted by the closure relations for a cooling 
break (jZhang Meszarosll2004l ). The Galactic X-ray absorption and reddening we re fixed to the 



values corresponding to Nh=1.0x 10 cm and E(B — U)=0.18 (lOates et al.ll201ll ). respectively. 



The intrinsic X-ray absorption was modeled by assuming an absorber with solar metallicity. To 
model the host intrinsic extinction we tested each of the Milky Way (MW), Lar ge Magella nic Cloud 
(LMC) or Small Magellanic Cloud (SMC) extinction law, as parametrized by iPeil (|1992l ). Limited 
statistics of our data did not permit any of these laws to be excluded or preferred. 

For epochs I and II, the SEDs are well fit by a single power law spectrum over the whole energy 
range from the keV (BAT data) to the GeV (LAT data), as shown in Fig. [5j 

For epoch III the analysis was complicated by the presence of X-ray flares and also limited by 
the poor optical coverage at early times (c.f. Fig. (H). In order to estimate the afterglow underlying 
the observed flaring activity, we extracted the X-ray spectrum in interval [160 s, 227 s] and rescaled 
it by extrapolating the light curve best-fit model. We extracted the spectrum during the X-ray 
flares from 56 s to 150 s, and subtracted from the estimated afterglow contribution. 

We also explored whether the observed S EP evolution was con sistent with the expected after- 
glow behavior in a wind-like density profile (jGranot Sarill2002l ). To this aim, we performed a 
joint fit by tying some parameters between the different epochs: we constrained the break energy 
to increase with time as E^ oc t 0,5 , and also held fix the host absorption and extinction as they 
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100 1000 10 4 10 5 10 6 

Energy [keV] 

Fig. 5. — SEDs at epochs I and II. GBM BGO data upper limits at both epochs are very similar, 
so for clarity only the upper limit at epoch I is shown. The SEDs at epochs I and II are well fit 
by a power law with photon indices, 1.871q ?! an d 1.95^0^, respectively. The continuous lines 
represent the best-fit, and the dashed lines the 1 a ranges. 
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are not expected to vary. We assumed a realistic model for the afterglow spectral shape, that is a 
smoothly broken power law of the form: 



F„(E) oc (E/E bk y s ^ + (E/E 



bk ) 



-1/s 



(6) 



where the curvature parameter s was held fixed at 0.8. Epoch III was excluded from the fit 
procedure, but compared to the resulting best-fit model. SEDs and best-fits using the wind model 
for epochs III, IV and V are shown in Fig. El and the best-fit results reported in Tab. [3] for all five 
epochs. 



-i iiiu| — i 1 1 iiinj — i 1 1 iiinj — i 1 1 1 1 1 n| — i 1 1 iniij — i 1 1 1 1 1 n| — i 1 1 iniij — i 1 1 iniij — i 1 1 1 ini| — i 1 1 iniij — i 1 1 1 1 1 n| — i 1 1 1 in 



III: t=100 s 
550 s 
7d 




io- 3 0.01 0.1 



10 100 1000 10 4 10 5 10 6 
Energy [keV] 



Fig. 6. — SEDs for epochs III, IV and V. For epochs III and IV, the spectra are built using UVOT 
(~3x 1(T 3 keV), XRT (~ 0.5 - 10 keV) and LAT (10 5 -10 7 keV) data. For epoch V, we have used 
the late time GROND observations (~ 3 x 10~ 3 keV) and XRT data. The solid lines show the fit 
of the wind model discussed in the text at each time, for the parameters reported in Tab. [3l 
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Model/Epoch 



1(9 s) 



II (16 s) III (100 s) IV (550 s) V (2.7 d) 



power law 
Index 



-1.87 



-0.11 
-0.07 



-1 95 +an 



Broken power law 

Popt 

Ebreak (keV) 

s (curvature) 
X 2 /d.o.f 



-0.45^ 



hO.07 
-0.09 

-0 Q5+ - 07 

u - yo -0.09 

04 +0 - 03 

u - uq: -0.01 



-0.66 
-1.16 



0.03 
0.03 
0.03 
-0.03 



0.8 



0.8 



- 01 
0.8 



654/731 



Table 3: Results of broadband fits of the SEDs. For epochs I and II the simple power law is the 
best model. Epochs IV and V have been fitted together with a smoothly broken power law, with a 
curvature s = 0.8, Ebreak ~ i ' 5 and /3 opt =f3 x +0-5. Epoch III was not part of the fitting procedure 
but we check that the data were compatible with the best-fit model obtained for epochs IV and V. 



5. Discussion and Interpretation 

5.1. Prompt emission 

GRB 110731A was bright in the LAT data with its most prominent peak in interval c. The onset 
of the LAT emission is delayed by 2.5 seconds with respect to the GBM To, a time delay comparable 
to the few-seconds delay observed in most long-duration LAT GRBs. The time-resolved spectral 
analysis suggests that an additional, hard emission component becomes increasingly prominent 
with time. Indeed, in b, the BAND+PL model fits the spectrum better than a simple BAND, 
with a significance only slightly below our threshold. While an additional component cannot be 
statistically resolved in c, an additional power law component improves the BAND-only fit in d 
as well, although with a rather low significance. By adopting the BAND+PL model for intervals 
b and d we can measure the ratio of the fluences of the power-law and BAND components, both 
measured in the 10 keV - 10 GeV energy range. The ratio is 20% in interval b and 35% in interval 
d, respectively. This trend confirms an emerging picture, corroborated by spectral analysis of other 
bright LAT bursts wit h evidence of a pow er law compone nt in addition to th e BAND component 
(namely GRB 090510. kbdo et all bold : GRB 090902B, kbdo et al.1 boOQal and GRB 090926 A, 



Ackermann et al 



2011 



, for which the flux of the additional power law component grows with time, 
and most likely dominates in the temporally extended LAT emission beyond the prompt emission 
phase. 

The ori gin of this additional spectral component is not yet fully understood. An early after- 



glow model (IKumar fe B arniol Duran 



De Pasquale et al 



2010 



Corsi et al 



2009; 



20101 : iRazzaque 



2010; 



Ghisellini et al. 



2010 



Kumar fc Barniol Duran 

201ol ). that is consistent with the above sce- 



nario, produces the power law component from the forward s hock of the GRB blast wave that 



propagates into the external medium surrounding the GRB (jMeszaros h Reesl 119971 ; ISari et al 
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19981 ). The delayed onset of the LAT emission is explained as the time required for the for- 
ward shock emission to become detectable in this scenario. In the context of the internal shock 



model (|Rees fc Meszaroslll994l ). the additional spectral component can ar i se du e to Compton scat- 
tering of soft targe t photons by rela tivistic electrons (jMeszaros Sz Reed 11994 ; IWang et al.1 120091 ; 
Bosniak et alj 120091 : iToma et al.l 1201 ll ). The 2.44 s delay for the LAT-detected emission in these 
scenarios would indicate GeV emission to be variable on a similar time scale and arising due to 
Compton emission from late internal shocks. Finally, hadronic emission, either proton/ion syn- 
chrotron ra diation or photo-pion-induced cascade radiation, can pro duce an additional spectral 
component (jAsano et al.ll2009i ; iRazzaque et al.ll2010l ; IWang et al.ll2009i ). 



The delayed onset of the high energy emission in the hadronic models is interpreted as due 
to the time required for proton/ion acceleration and cooling, as well as to the time required to 
form cascades. In particular, to generate the 2.44 s delay in LAT the jet magnetic field needs to 
be B' « 10 5 (r /500)- 1 / 3 (t onset /2.44 s)- 2 / 3 (£ 7 /100 MeV)" 1 / 3 G, with a jet bulk Lorentz factor 
To, in case of proton-synchrotron model and the corresponding isotropic-e quivalent jet luminos- 



ity is L iBt > 10 57 (r n /500) 16 / 3 (t nTlsP ,/2.44 s) 2 / 3 (£ 7 /100 MeV)" 



2/3 er g g i (jRazzaque et al 



2010 



Wang et al.l [2009J) . While the hadronic models require a larger total energy than the leptonic 
models ; the energy budget can be brought dow n to acceptable levels for a bulk Lorentz factor 

In addition, the flat spectrum (index ~ —2) for 



< 500 ( Abdo et al 



2010 



Ackermann et al 



2011 



the extra component is also favorable to the energy budget, and the observed low-energy (10 keV) 
excess (c.f. Fig. [3j top) could be additional evidence for a hadronic model. These tw o properties 



are very similar to the f eatures already reported for the LAT burst GRB 090902B (lAbdo et al 
2009al : ksano et al.lhoioh . 



At a redshift z = 2.83, the peak isotropic-equivalent luminosity (10 keV-10 GeV) of GRB 110731A, 
as measured in interval c, is A so , p k = 6-3i J X 10 53 erg s _1 . The corresponding peak of the energy 
spectrum is at E^ = (2 + T a )Eo = 683^ 2 gQ keV, from the BAND fit. The total isotropic-equivalent 
energy (10 keV-10 GeV) over the full prompt phase (PI) is Ei so = 6.01q 2 x 10 53 erg, assuming 
a COMP+PL model. The BAND+PL fit to the same spectrum results in a somewhat larger 
E iso = 7.6±g;| x 10 53 erg with a BAN D £ pk = (2 + rjfi ~ 300 keV. These values are broadly con- 
sistent with the empirical £iso,pk~-Epk (jYonetoku et al.ll2004l ) and Ei so -E p k (lAmati et al.ll2009i ) rela- 



tions for long GRBs. Another empirical relation between the luminosity and the spec tral lags that 
we measured to be —41 ±28 ms between two Swift BAT energy bands, is also satisfied (jNorris et al 



20001 : lUkwatta et aUl2012 ) 



Constraining the bulk Lorentz factor, Tq, of the jet is a major challenge in GRB science. A 
constraint from the 77 pair production opacity argument can be used to derive a lo wer limit on Tp 
if the spectrum does not have a cutoff and extends t o the highest observed energy (|Krolik k, Pier 
199ll : iFenimore etaDll993l : baring fc Harding! Il997l : iLithwick fc SarilboOll : IRazzaque etall 120041 ) . 
For GRB 110731A, a 2.0 GeV photon was detected at 8.27 s. The spectrum in interval d which 
included this photon, is preferably fit with a BAND+PL. We derive a lower limit, r m ; n ~ 600, 
from r 77 < 1 using the above information and the median value of the FWHM of the pulses within 
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the IgOjGBM; At = 0.35 ± 0.02 s, as the variability time in the single-zone emission ap proximation 



with a homogeneous distribution of photons (for a description of the method, see lAbdo et al 



201u ) . The corresponding emission radius is R > r^ in cAt(l + z)' 1 ~ 10 15 cm. In interval P2 we 



detect a high-energy cutoff with a significance only slightly below our threshold. This provides 
an opportun ity to estimate Tn as in GRB 090926A, by assuming that the cutoff is due to 77 pair 



production (jAckermann et al.ll201ll ). We derive T = 530 ± 10, using the COMP+COMP model 
with a 390^20 MeV folding energy for the high-energy component, and with At = 0.43 ±0.03 s, the 
median value of the FWHM of the pulses in P2. The maximum photon energy in this time interval 
is 0.9 GeV. Observe that while regions somewhat smaller than those inferred from these choices of 
At can be possible, since the pair opacity is a strong function of Tq, the values of the bulk Lorentz 
factor inferred here are only weakly depen dent on the choice of At. Note that effects within the 
emission region such as radiation transport (lAbdo et al.ll201Cf) o r a time-dependent increase of r 77 



and geometrical effects (jGranot et al 



2008; 



Hascoet et al 



2012j) can reduce the value of r m i n or Tq 
by a factor of up to 3, even for a single-zone mo del. Indeed the for the interval d reduces to a 
lower limit of ~ 300 for the parametrization by lHascoet et al.1 (|2012l ). Limits on the bulk Lorentz 
factor can be further relaxed for two-z one emission models, where GeV ph otons are emitted from 
a larger radius than the MeV photons (IZou et al.ll201ll : lHascoet et al.ll2012l ). Such reductions in Tq 
clearly can prove important for the viability of hadronic models, as discussed above. 



5.2. LAT temporally extended emission and mult iwave length afterglow 



Thanks to the wealth of simultaneous multiwavelength data, we can test various afterglow 
models and extract the preferred model parameters. In particular, these observations can help 
constrain whether the temporally extended LAT emission originates in the forward shock that 
produces the late-time afterglow emission. 

We discuss epoch IV first (c.f. Fig. Hj) since the data are most constraining in this time 
interval. In Section [3.21 we observed that the optical flux (a op t = 1-37 ± 0.03) decays faster than 
the X-ray flux (a x ,i = 1.10 ± 0.02 or a x = 1-189 ± 0.007), and in Section we reported 

that the broadband SED is well fit by a smoothly broken power law at epo ch IV. Both features 
favor a wind rather than I SM afterglow model, with a slow-cooling spectrum ([Chevalier &: Lill2000l : 
Panaitescu &: Kumarll2000l ). This results from an examination of the relations between the decay 
and spectral indices t hat describe flux evolution, F v oc t~ a v~^ , at different energies of the afterglow 



synchrotron spectra (jGranot &: Sari 



2002 



Zhang et alJ l2006). The break in the SED at epoch IV 



can be readily interpreted as a cooling break (u c ) rather than a break due to the typical synchrotron 
frequency of the minimum energy electrons (f m ), in both the slow- and fast-cooling spectra; the 
respective synchrotron indices below such breaks are /3 = 2/3 (slow-cooling) or /3 = 3/2 (fast- 
cooling). This is because the fast decay of the observed optical flux is contrary to (i) an increasing 
(ISM) or invariant (wind) flux behavior expected below v m in the slow-cooling spectrum; and (ii) 
a much slower decay behavior, F v oc t -1 / 4 (ISM or wind), expected below v m in the fast-cooling 
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spectrum. Further investigation shows that the fast-cooling spectrum is disfavored at epoch IV, 
since (i) ax,i or ax and px,i (= 0.87±0.04, derived as Tx — 1 before the temporal break) for the X- 
ray flux are both incompatible with the F v oc i -1 / 4 ^ -1 / 2 behavior expected for v c < vx < v m (ISM 
or wind); and (ii) a op t is incompatible with either F v oc i 1 / 6 (ISM) or F v oc i~ 2 / 3 (wind) behavior 
for u opt < v c . The spectral index /3x,sed = 0.95^q'q7 for v c < ux, fitting the broadband SED at 
epoch IV, is compatible, within 1-2 a, with f3x,i and with the expected /3 = (2/3)a + 1/3 behavior 
in the slow-cooling spectrum for ax,i or ax- While the above is common behavior for both the ISM 
and wind models, the observed a op t favors a (3 = (2/3)a — 1/3 relation (wind) over a (3 = (2/3)a 
relation (ISM) when compared with the broadband SED fit result /3 pt,SED = /3y,sed — 0.5. The 
expected optical spectral index is within 2 a of /3 op t,sED in the wind model and deviates by more 
than 5 a from /3 pt,SED in the ISM model. Thus our analysis favors a slow-cooling spectrum and 
wind environment for broadband data at epoch IV. 

The broadband SED at 2.74 d, epoch V, has far less predictive power than the SED at epoch 
IV, with a single power law being statistically favored over a broken power law. Nevertheless, 
the joint fit to both SEDs using a smoothly broken power law, with tied indices but allowing 
the break frequency to increase as oc t 1 / 2 between the SEDs as expected in the wind model, re- 
sulted in a good fit. This result further supports the wind model. A steeper decrease in the 
X-ray flux (ax,2 = 1-32 ± 0.03) after 4.6l^'g ks can be explained as the break frequency v c ap- 
proaching the X-ray band with time and reaching ~ 1 keV at 2.74 d. An optical rebrightening 
is clearly visible in GROND data at epoch V w ithout an accompanying rebrightening in X-ray, 
as observed in GRB 081029 (jNardini et al.ll201ll ). The optical rebrighte ning could possibly orig- 
inate from a t wo-component jet st ructure as modeled in GRB 030329 (IBerger et al.1 l2003l ) and 
GRB 080319B (IRacusin et al.ll2008l ). However, the lack of data during the rise of the optical flux 
and sparse late-time data prohibit further study of the optical rebrightening in GRB 110731 A. 

The LAT temporally extended emission that is contemporaneous with the BAT and GBM 
detected emission at very early times, epochs I and II, can be fit with single power-law models to 
produce a broadband SED. The resulting photon indices, 1 .87^0'^ and 1.95+°;^ respectively, are 
compatible with the measured X-ray photon index Tx at later times. In the context of the wind 
afterglow model, these two very early SEDs can be interpreted as being due to emission above 
v m in a fast-cooling spectrum [y c < v m ) which subsequently turns into a slow-cooling spectrum 
{y m < v c ) as v c increases with time and moves past v m , giving rise to the broken power law SED 
at epoch IV. However, while the LAT flux decay index «lat = 1-55 ±0.20 is marginally consistent 
within 2 a with the later X-ray flux decay index ax,i or ax, the BAT flux decays (obat = 2.3±0.3) 
at a much faster rate. A steeper decay of the BAT flux can be due to an emission component, suc h 
as high-latitude (9 > Tg 1 ) emission from the fireball ( Kumar &: Panaitescul 120001 ; iDermerl 12004 ) . 



in addition to the underlying afterglow emission. Indeed, fitting the BAT flux as F v oc (t — i p k)~ a , 
where t p k = 5.5 s was chosen to coincide with the brightest peak in the LAT light curve, results in a 
flux decay index (a;BAT,t pk = 1-25 ± 0.15) which is more compatible with «lat and ax,i or ax, yet 
captures the steep decay of the BAT flux at very early time. Part of the LAT emission at ~5.5 s 
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likely originates from such high-latitude emission and/or internal shocks. Indeed an extrapolation 
of the LAT afterg low flux at this time fails to reproduce all of the observed emission, as also noted 
for GRB 090510 fae et al.l boilf ) . 



Epoch III is rather complex with 2 weak X-ray flares. The most likely origin of the flares and 
a rebrightening of the emission detected by the B AT in coincidence with the first flare is late-time 
activity of the central engine (jZhang et al.l 120061 ) , as also implied by the noticeable variation of 
the spectrum from the underlying afterglow emission. The change in the afterglow synchrotron 
spectrum from the fast- to slow-cooling regime is expected to take place in this interval or earlier, 
since the SED in the next interval is best fitted by a broken power law as mentioned above. 

Following the above discussion on the broadband SED of GRB 110731A for multiple epochs 
and flux decay behavior, two main features can be highlighted: 

• LAT temporally extended emission, as early as 8.3 s as seen in the first broadband SED, is 
compatible with the afterglow synchrotron emission in other bands. 

• A wind afterglow model is favored over an ISM model, from the behavior of the SEDs. 



In the wind afterglow scenario, however, the onset of the afterglow needs to take place quite rapidly. 
There are several hints that this is the case: (i) the brightest peak in the LAT light curve dominating 
the vF v flux at all frequencies takes place as early as 5.5 s; (ii) the absence of a brighter flux in the 
BAT or GBM at a later time; and (iii) the LAT flux decays smoothly after 5.5 s. While the last 
point is quite obvious, the first two hints result from the fact that th e bolometric flux of the forward- 
shock emission peaks at the deceleration time of the GRB fireball (IBlandford &: McKedll976l : ISari 
19971 : iGhisellini et al.ll2010l ). In the case of an adiabatic fireball in the wind environment, the 



deceleration time is idee = (1 + z)Ek/ '(16TTm p c 3 ATq), where Ek is the isotropic-equivalent kinetic 
energy of the fireball during the decelerating phase and A = 3.02 x 10 35 j4* cm -1 is the wind 
parameter for a 10~ 5 M(n yr -1 mass-loss ra te in the wind of velocity 10 3 km s , with A* = 1 as 
an scaling parameter ([Chevalier Lill2000i ). Note that formally the afterglow onset requires that 
^dec ^ ^GRBj where Tqrb is the duration of the prompt phase, for the "thin-shell" formula used 
here. Nevertheless, most of the emission detected by the GBM (50-300 keV) comes within the 
^95,GBM (7.6 ± 0.3 s), which is similar to the deceleration time scale needed to explain the earliest 
broadband SEDs. Moreover, a fraction of the emission arriving later in the formal prompt phase 
may originate from the very early afterglow. 



5.3. Afterglow parameters in the wind model 

Adopting the wind afterglow model and the broadband SED fit parameters of epoch IV, 
we derive the afterglow model parameters for GRB 110731 A. The optical and X-ray flux decay 
indices constrain the electron acceleration index to be p = (4/3)a op t + 1/3 = 2.16 =L 0.04 and 



- 30 - 



p = (4/3)ax,i + 2/3 = 2.13 ± 0.03, respectively, from the F v oc t~ Q V _/3 relations. For the sak e of 
simplification, we assume p = 2.2. We then use the parametrization by iGranot Saril (|2002j) for 
the break frequencies (u m , v c ), the fluxes at the break frequencies (F Um , F Uc ), and the transition 
time (to) from the fast- to slow-cooling spectrum as 



hv m 
hv c 

to 



3.5 x 10 4 e 2 M 2 E^tf 12 keV, 



2.8 
7.1 
1.2 



x 10- 8 e^A^EUS 1 / 2 keV, 



10 



.1/2 



ln 12 1.2 1.7 42.2 pl/2. -1.7 
1U e e £ B A^ £/ 55 



mJy, 



1.1 x 10 7 e e e B A* s. 



(7) 



Here e e and e b are the usual micro-physical parameters (ISari et al.1 Il998l : IGranot k, Saril 120021 ) , 

^55 = E k /W 55 erg, and t\ = i/10 s. 

Using the reference time t c ~ 550 s for the cooling break in the SED, hu c ~ 0.04 keV and 
F Vc ~ 0.5 mJy, we obtain acceptable afterglow parameters for i^ec ~ 5.9 s, the middle of the time 
interval with the brightest LAT peak, as: 



E k = 3.3 x 10 



51 



3.2 x 10" 



A. 

0.05 



-5/6 



tde 



5.9 s 



-8/9 



500 J 



-32/9 



tr 



17/18 



e B = 1.1 x 10 



-2 



0.05 



550 s 

5.9 sy Vsooy 



A* 
0.05 



0.04 keV 



5.9 s 

17/18 



5lV 



erg, 



4/3 



550 s 



to = 19.5 



0.05 



-5/6 



tde 



5.9 s 



-5/9 



(IL 
\500 



-20/9 



t, 



550 s 



23/18 



hv. 



0.04 keV 



1/3 



5/18 



500 J 
0.5 mJy / 



0.04 keV 



0.5 mJy 



5/6 



s. (8) 



Thus the fast- to slow-cooling transition takes place towards the end of epoch II, in which a single 
power law is a better fit to the broadband SED. The coasting bulk Lorentz factor T$ ~ 500 is 
compatible with the value/lower limit obtained from the 77 opacity argument for the internal 
shocks. The kinetic energy is also a factor ~ 5 larger than the isotropic-equivalent 7-ray energy 
Eiso = (6.8 ± 0.1) x 10 53 erg, which is compatible with the typically assumed ratio between the two 
energies. A deceleration time later than 5.9 s but earlier than ~55 s can also be accommodated 
with acceptable parameter values in this scenario, without large deviations from the Tq ~ 500 and 
A* ~ 0.05 values used in Eq. ([8]). However, the fact that the broadband SED at epochs I and II 
can be fitted with power laws, strongly argues in favor of an afterglow onset time before 8.3 s. The 
variations of the parameter values derived in Eq. are minimal in this case. With the parameters 
E k and A+ fixed, a change of td ec from 5.9 s to 8.3 s corresponds to Tq ~ 545. 



The maximum energy of the synchrotron photons, which is independent of the micro-physical 
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parameters (see, e.g.. iRazzaque et al.ll2010l ). can be written as: 
/^ max « 0.24</»- 1 (l + z)- 1 r GeV 

A \ -1/4 

15 



-U-l/4 

r 2 



0.05 



St 



10 54.5 



erg 



1/4 



GeV. 



(9) 



Here < 1 is the acceleration efficiency for electrons and £2 = t/100 s. This equation applies 
to scenarios where lepton acceleration is gyroresonant and is radiation-reaction limited by syn- 
chrotron cooling. In such cases, the maximum electron energy in the comoving jet frame is of the 
order of m fi c 2 divided by the fine structure constant, a result well-known i n the context of active 
galaxies (jGuilbert et al.lll983l ) and the Crab nebula ()De Jager et al.lll996l ). Thus detection of a 
3.4 GeV photon at ~436 s from GRB 110731 A (see end of Section 14. 1 . 1 [) is consistent with the 
maximum synchrotron photon energy for the afterglow parameters derived in Eq. ([8]). 



5.4. Comparisons with other GRBs 



The sample of LAT-detected GRBs is still relatively small, with a detection rate of ~ 8/yr (jPiron 



20121 ). Here we briefly comment on the properties of GRB 110731A in comparison with other LAT- 
detected GRBs that were bright at GeV energies, and with GBM-detected GRBs in general. More 
detailed comparisons will be presented in a catalog of LAT GRBs that is in preparation. We also 
comment on the X-ray and optical afterglow of GRB 110731 A in comparison with other GRBs. 

The total (10 keV-10 GeV) isotropic-equivalent energy i£ 7( i SO of GRB 110731 A falls in the 
middle of the distribution for LAT bursts with known redshifts; only GRB 100414A, GRB 091003, 
GRB 090328 and the short burst GRB 090510 have lower £ 7jiso . In the 1 keV-10 MeV range, E J>iso 
for GRB 110731A is lower than for a number of GBM-only GRBs in the same redshift range. The 
fact that GRB 110731A was very close to the LAT boresight played a crucial role for its detection 
at > 100 MeV. An additional power law component as in GRB 110731 A has been detected in all 
bright LAT-detected GRBs and a cutoff in the power law for the > 100 MeV energy range has 



been reported in GRB 090926A (jAckermann et al.ll201ll ) in addition to the present case. The value 



for simple single-zone model, in the range of > 200 to > 1000 for bright LAT GRBs ( 


Abdo et al. 


2009c 


Greiner et al. 


2009; 


Abdo et al. 


2009a 


2010: 


Acker mann et al. 


2011 


Cenko et al. 


2011 


). The 



To ~ 500 value that we derived for GRB 110731 A using the same argument is in the middle of this 
range. Note that uncertainties in underlying emission modeling may scale down these values by a 
factor 2-3 (Sec. 15. ip . Although the current sample is rather small, a cutoff in the BAND or PL 
spectrum could be one reason for the LAT non -detection of a number of bright GBM bursts within 
the LAT field of view (jAckermann et alJboia ). 



Temporally extended high-energy emission detected in GRB 110731 A is quite common among 
LAT-detected GRBs, showing a smooth power law decay of the flux. T he flux decay index, ho wever, 



is not the same among bursts but varies between ~ 1 and ~ 1.5 (jGhisellini et al.1 l2010l ). The 
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common behavior of the temporally extended LAT flux is compatible with afterglow emission, but 
a direct model comparison w ith contemporaneous multiwavelength dat a sets is possible only for 
the short burst GRB 090510 (|De Pasquale et alJboiOl : Ubdo et alJboiOl ) and for the current long 
burst GRB 110731A. 

The X-ray afterglow of GRB 1 10731A does not display the "canonical" steep-, flat-, normal- 
decay behavior (jNousek et al.ll2006l ) observed in ~40% of all long GRBs detected by Swift. This 
trend for GRB 110731 A is also o bserved in other L AT bursts with earl y-time Swift XR T ob- 
servations such as GRB 100728A (|Abdo et alJboilh and GRB 110625A (|Tam et alJ boid l. and 
may suggest that the afterglows of LAT bursts tend to be dominated by the bright forward-shock 
emission rather than prolonged episodes of energy injection. 



We have studied the optical afterglow of GRB 1107 31 A using the method of Kami et alJ (12006) 
and compared it with other well-measured afterglows ( Kann et al. 2006 ; Kann et alJ 201ol . 2011 ). 
As the data quality in the optical is not sufficient to allow a fit with all parameters free, we fix the 
underlying spectral slope to the value derived from the broadband X-ray-to-optical fit, = 0.66. 
A free fit to the optical data alone results in a similar value, ~ 0.7, but with very large errors. 
Assuming an SMC-like dust density, we derive an extinction Ay = 0.24 ± 0.06 (note that this 
error is somewhat underestimated due to the constrained spectral slope). Alternatively, without 
dust, the optical data yield a red spectral slo pe of = 1.41 ± .18. Using this extinction value and 
the redshift, we find a magnitude correction ( Kann et al.ll2006l ) of dRc = — 3.46i ' 2 3 (alternatively, 
dRc = —3.09 ±0.13 without extinction). Using the known SEP, we are a ble to create a com pound 

The 



2011 



Tanvir et al 



2011 



/2-band light curve, to which we add GCN data (IMalesani et al. 
resulting light curve spans from 32 seconds to « 1.5 days after the GRB, and reaches a peak 
at 9 th magnitude, implying that this is initially an exceptionally bright afterglow. Comparing 
the known sample of optical afterglows with prompt-emis sion detections in the high-MeV/GeV 



range (including one AGILE GRID GR B, GRB 080514B, iGiuliani et alJ 120081 : iRossi et al 



McBreen et al. 



2008 



2010 



Cenko et al 



2011 



we find that while some afterglows are among the most 
luminous known (GRB 090923, GRB 090926A), the others are of average luminosity at t = 1 day, 
including GRB 110731 A. The span is four magnitudes, an d if one extrapolates the steep late decay 
of the short GRB 090510 (jNicuesa Guelbenzu et al.ll2012l ). the span reaches about ten magnitudes. 
The very early detection of the afterglow of GRB 110731 A also allows us for the first time to 
compare the brightness at 0.001 days (normalized at z = 1) for a LAT-detected GRB, we find 
R = 10.5. Coupled with the unbroken decay from the ea rliest detection on , this makes it one of 
the brightest forward-shock dominated afterglows known (|Kann et al.ll2010l ). 



Our broadband spectral and temporal modeling, including LAT data, favors a circum-burst 
medium with a wind-like density profile and an afterglow onset time < 9 s after trigger. This is 
the first time that a LAT-detected long GRB could have such an early onset time estimated using 
multiwavelength data. The corresponding bulk Lorentz factor value T$ ~ 500, deduced from the 
deceleration time scale of the fireball, is also better motivated for GRB 110731 A. With a constant 
density ISM environment and by assuming that the deceleration time scale coincides with the peak 
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LAT flux time at 0.7 s, the initial bulk Lq rentz factor of the short GRB 090510 was calculated 
to be To ~ 2000 (|De Pasquale et al.l |2010| ) . Using Tq above the lower li mit T m m, calculated 
from 77-opacity argument in the simple one-zone model in the prompt phase, ICenko et al.l (|201ll ) 
have modeled the temporally extended X-ray-optical-radio afterglow data of 4 LAT-detected long 
GRBs: GRB 090323, GRB 090328, GRB 090902B and GRB 090926A. A wind environment was 
preferred for each of these G RBs except for GRB 090902B, for which an ISM environment was 
preferred (jCenko et al.l l201ll ). Although the sample is rather small, it hints at an interesting 



trend since a systematic study of well-sampled long G RB afterglows show s that the majority of 

The micro-physical 



them is consistent with a uniform density environment (jSchulze et al.ll201ll ) 
parameters, e e and €b that we estimated for GRB 110731A in Section 15.31 are lower than for 
GRB 090510 with e e ~ €b ~ .1 which were also estim ated from broadband spectral and temporal 
modeling including L AT data (|De Pasquale et al.ll2010l ) . Value s of e e and en close r to our estimates 



have been obtained for GRB 090323 with late afterglow data (jCenko et al. 



2011 



6. Conclusions 

GRB 110731 A is the first long burst detected by instruments on-board Fermi and Swift with 
simultaneous coverage from optical to gamma rays for a few hundred seconds. Prompt follow up 
observations by MOA and late-time data from GROND provide crucial optical data from a few 
minutes to a few days after the trigger and while the afterglow is still bright in the XRT. The 
redshift of the burst z = 2.83 is relatively high but within the range of other bright long bursts 
detected by Fermi LAT. 

GRB 110731A has a rather high total energy output, with an isotropic-equivalent gamma-ray 
(10 keV-10 GeV) energy release of E- lso ~ 6 x 10 53 erg within the prompt phase of the first ~ 8.5 s. 
The brightest peak in the LAT light curve at 5.5 s coincides with peaks in all other energy bands, 
producing a peak isotropic-equivalent 7-ray (10 keV-10 GeV) luminosity of ~ 6 x 10 53 erg s _1 . The 
LAT emission is delayed by ~ 2.5 s from the GBM trigger and continues for the next ~ 850 s. These 
features (delayed onset and temporally extended emission in the LAT) are common to many bright 
LAT bursts. Similarly to other bright LAT bursts, we detect an additional power law component to 
the Band function in the time-integrated spectral analysis, with evidence also in the time-resolved 
analysis, although with lower significance. For GRB 110731A, we also have evidence for a cutoff 
in the power law component at GeV energies, although with a significance which is not as high 
as that measured for GRB 090926A. Indeed, two Comptonized models provide a slightly better fit 
to the time-integrated spectra from Tj^lat to Tgs^BM with respect to the Band plus power law 
model. Using a measured variability time scale in the GBM data and by assuming that the LAT 
photons are co-spatial with softer photons, we calculate a jet bulk Lorentz factor of ~ 500 from 
the 77 attenuation mechanism producing the GeV spectral cutoff. A value smaller by a factor 
~ 2 to 3 is also acceptable within modeling uncertainties; for instance, we found ~ 300 using a 
parametrization for two-zone emission models (c.f. Sec. 15 . X [) . 
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The broadband spectrum of GRB 110731A from optical to > 100 MeV at 550 s, when fit with 
a broken power law, together with temporal flux decay behavior in different wavelengths, favors an 
afterglow wind model as the origin. We interpret the spectral steepening as a cooling break and 
verify its expected temporal evolution by fitting the broadband spectrum at 2.74 d. The presence 
of two mild X-ray flares and lack of good optical data at an earlier epoch (III) prohibit us from 
independently verifying the afterglow model. However, the broadband emission, after extraction of 
the X-ray flares, is consistent with the model at 550 s. Most remarkably, we find that BAT and 
LAT spectra at 9 s and 16 s can be fit with single power laws with compatible spectral indices, which 
in turn are consistent with the measured X-ray spectral index at later times. These broadband 
spectra at early times strongly suggest their origin as afterglow emission as well. This scenario, 
however, requires an afterglow onset time before ~ 8 s. 

With an afterglow onset time of 5.9 s, first time interval after the brightest LAT peak, we 
interpret subsequent multiwavelength data as originating from forward shock emission in a wind- 
type medium. We cannot rule out contamination starting at an early time from prompt emission, 
however. With an initial bulk Lorentz factor ~ 500-550 and wind density parameter within the 
acceptable range we are able to fit and interpret the broadband data. It is worth emphasizing that 
this second derivation of the bulk Lorentz factor, using the afterglow modeling, is independent from 
the estimate done using the cutoff during the prompt phase, but that both numbers are nevertheless 
very compatible. The resulting micro-physical parameters are e e ~ 3 x 10~ 3 and e# ~ 10 -2 which 
favor a somewhat larger magnetic energy density than electron energy density in the forward shock. 
The time scale for the fast-to-slow cooling spectral transition is also compatible with earlier single 
power law spectra and later broken power law spectra. The total isotropic-equivalent jet kinetic 
energy calculated is ~ 3 x 10 54 erg which puts GRB 110731 A among the most energetic bursts. 

Compared to other LAT GRBs, GRB 110731 A is in the middle of the fluence distribution and 
is not as bright as the four brightest GRBs, but due to favorable observing conditions we were 
able to collect comprehensive multiwavelength data and put strong constraints on the temporally 
extended GeV emission. In addition, both X-ray and optical data strongly suggest the presence of 
early forward-shock emission in this burst, which would constitute the earliest afterglow emission 
detected so far among LAT-detected GRBs. 

We find that temporally extended LAT emission is compatible with originating from the for- 
ward shock that produces the broadband afterglow emission. In the context of this scenario, the 
delayed onset of the LAT emission can be interpreted as the increasing flux of the forward-shock 
component, before reaching the deceleration time and after which the traditional afterglow phase 
begins. The validity of this scenario and the origin of the LAT emission will be further tested by 
future multiwavelength observations of Fermi LAT-detected bursts. 
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A. 


Data tables 




Time 


Energy 


Photon Flux above 100 MeV 


Test Statistic 


Bins (s) 


Index 


(ph cm _2 s _1 ) 




(11 2 35-3 59 


-2.22 ± 0.39 


1.11 ± 0.40 x 10~ 3 


73 


(2) 3.59-4.56 


-2.83 ± 0.60 


1.47 ± 0.55 x 10~ 3 


53 


(3) 4.56-5.47 


-2.45 ± 0.49 


1.45 ± 0.57 x 10~ 3 


68 


(4) 5.47-5.67 


-2.25 ± 0.41 


6.80 ± 2.48 x 10~ 3 


88 


(5) 5.67-6.58 


-3.12 ± 0.68 


1.69 ± 0.60 x 10~ 3 


90 


(6) 6.58-8.27 


-2.11 ±0.37 


7.66 ± 2.87 x 10" 4 


70 


(7) 8.27-11.54 


-1.80 ±0.31 


3.50 ± 1.42 x 10~ 4 


73 


(8) 11.54-19.73 


-4.65 ± 1.24 


1.92 ±0.71 x 10~ 4 


41 


(9) 19.73-32.60 


-3.13 ±0.82 


8.39 ± 3.68 x 10~ 5 


27 


(10) 32.60-110.97 


-1.92 ±0.38 


1.24 ±0.65 x 10~ 5 


20 


(11) 110.97-227.04 


-3.01 ±0.68 


1.64 ±0.71 x 10~ 5 


17 


(12) 227.04-853.89 


-1.69 ±0.35 


1.77 ± 1.16 x 10" 6 


16 


(13) 853.89-1433.65 


-2.25 


< 3.9 x 10~ 6 


2 



Table 4: LAT Time-resolved spectroscopy data, photon fluxes and photon indices. The last entry 
is an upper-limit assuming a power-law index of —2.25. 
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Time (s) Exposure(s) Magnitude Flux (mjy) Filter 



62 
627 
800 


9 

20 
20 


13.081";* 

16.37+^22 
16.86l° : « 


21.26 ±0.71 
1.03 ±0.14 
0.66 ±0.13 


V 
V 
V 


1100 
12753 


619 

436 


17.481°;^ 
20.301^ 


0.37 ±0.06 
0.02 ±0.06 


V 
V 


1557278 


1069878 


> 20.42 


< 0.03 


V 


552 
725 


20 
20 


17-62l°-2° 
18.021°;™ 


0.36 ±0.06 
0.25 ±0.05 


b 
b 


1239 


192 


18.91+n'5? 
^ J — 0.41 


0.11 ±0.03 


b 


15599 


17857 


zi.io_ 73 


0.01 ±0.01 


b 


1557674 


1070041 


> 21.00 


< 0.02 


b 


299 


25 


17.28l°f 


0.18 ±0.02 


u 


324 
349 


25 
25 


17.051°,;! 

16.9818:1? 


0.22 ±0.02 
0.23 ±0.02 


u 
u 


374 


25 


16.7618$ 


0.28 ±0.02 


u 


399 
424 
449 


25 
25 
25 


17 n +0.29 

1 '- 11 -0.23 

17 co+0.39 
i(.oo_ 02g 

17 ^7+°- 41 


0.21 ±0.02 
0.14 ±0.02 
0.13 ±0.02 


u 
u 
u 


474 
499 


25 
25 


17 5Q+0.41 
17 ocr+0.34 


0.13 ±0.02 
0.16 ±0.02 


u 
u 


524 


25 


18.081°;^ 


0.08 ±0.02 


u 


1001 


620 


> 18.23 


0.01 ±0.01 


u 


21275 


29621 


> 20.37 


< 0.01 


u 


439746 


415954 


> 21.1 


< 0.01 


u 


976 


620 


> 18.14 


< 0.05 


uvwl 


951 


620 


> 19.12 


< 0.02 


uvm2 


989 


793 


> 18.1 


< 0.04 


uvw2 



Table 5: Swift UVOT data table for the individual filters. Time is the mid-time in exposure, in 
seconds, since BAT trigger. Exposure is the half-width of the integration duration in seconds. 
Magnitudes and flux densities have been corrected for Galactic extinction only. Upper limits are 
given at 3 a for both magnitudes and fluxes. 
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Time (s) Expo(s) Magnitude Flux (mjy) Filter 



80 


10 


- . . — 4-0 (17 


3.13 ±0.08 


white 


90 


10 


14.55i°l 7 


2.91 ±0.08 


white 


100 


10 


14.681°;°; 


2.59 ±0.07 


white 


110 


10 


14.87l°,g 


2.17 ±0.06 


white 


120 


10 


15.051q;q* 


1.83 ± 0.06 


white 


130 


10 


15.081°;°? 


1.79 ±0.05 


white 


140 


10 


1 'S 1 1+0.08 


1.74 ±0.05 


white 


150 


10 


1 r 9Q+ 08 


1.48 ±0.05 


white 


160 


10 


1 5 S6+ 008 


1.38 ±0.05 


white 


170 


10 


i r t;9+0.09 


1.19 ±0.04 


white 


180 


10 




1.22 ±0.04 


white 


190 


10 


1 5 6+°- 09 


1.11 ±0.04 


white 


200 


10 




0.98 ± 0.04 


white 


210 


10 


1 5 83+ ' 10 


0.89 ± 0.03 


white 


220 


10 


15 Q1+ 010 
io.yi_ 09 


0.83 ± 0.03 


white 


577 


20 


1 7 12+0.16 


0.27 ±0.02 


white 


750 


20 


17.53l°;?| 


0.19 ±0.02 


white 


940 


150 


18.2618:1| 


0.10 ±0.01 


white 


1178 


20 


19.32lJ|3 


0.04 ±0.01 


white 


1350 


20 


18.48l°;| 


0.08 ±0.01 


white 


6975 


200 


20.27l°;g 


0.01 ±0.01 


white 


104187 


14049 


> 20.96 


< 0.01 


white 



Table 6: Swift UVOT data table for the co-added exposures. Time is the mid-time in exposure, 
in seconds, since GBM trigger. Exposure is the half-width of the integration duration in seconds. 
Magnitudes and flux densities have been corrected for Galactic extinction only. Upper limits arc 
given at 3 a for both magnitudes and fluxes. 
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Time (s) Magnitude f Flux (mjy)t Filter 



Time (s) Magnitude f Flux (mjy)t Filter 



205 
268 
333 
428 
556 
656 
720 
785 
868 
933 
996 
1082 
1145 
1210 
1293 
1358 
1421 
1506 
1569 
1634 
1717 
1783 
1846 
1931 
1994 



14.50 ± 0.10 

14.87 ± 0.10 
15.22 ± 0.10 
16.76 ± 0.10 
17.19 ± 0.10 

16.29 ± 0.11 
16.42 ± 0.12 

16.54 ± 0.12 

17.88 ± 0.12 
17.97 ± 0.12 
17.99 ± 0.12 
17.09 ± 0.13 
17.16 ± 0.14 

17.30 ± 0.14 

18.52 ± 0.14 

18.55 ± 0.15 

18.53 ± 0.15 
17.68 ± 0.16 
17.81 ± 0.17 
17.70 ± 0.16 
18.84 ± 0.17 
18.97 ± 0.18 

18.89 ± 0.17 
18.03 ± 0.18 
18.03 ± 0.19 



4.1603 
2.960^ 
2.160j 
0.746j 
0.503j 



fO.408 
-0.372 
f 0.293 
-0.266 
f0.218 
-0.198 
f0.071 
J -0.065 
f0.051 
J -0.046 

O.8OO_ 080 
u -' ua -0.072 

o fisq+ - 074 
266+ - 031 

U.ZDD_q 028 

n 944+0-029 
U.^44_ Q 02 g 

-,+0.029 
-0.026 



0.2410 

0.384™ 
0.359 
0.315 
0.148 
0.143 
0.146 
0.222 



0.049 
-0.043 
0.041 
0.039 
0.021 
0.018 
0.021 
-0.018 
0.021 
-0.018 
0.035 
-0.031 



n 1 QS+0-033 
U ' ±JO -0.029 

0.219 
0.110 



0.098 
0.105 
0.162 
0.162 



0.03.-) 
0.030 
0.019 
0.016 
0.018 
0.015 
0.018 
-0.016 
0.030 
-0.025 
0.031 
-0.026 



I 
I 
I 

V 
V 

I 
I 
I 

V 
V 
V 

I 
I 
I 

V 
V 
V 

I 
I 
I 

V 
V 
V 

I 
I 



2060 
2143 
2208 
2271 
2362 
2426 
2489 
2573 
2636 
2787 
2851 
2914 
2999 
3127 
3214 
3277 
3342 
3637 
3978 
4061 
4126 
4188 
4485 
4548 



18.15 ± 0.19 

19.15 ± 0.20 
19.23 ± 0.20 
19.33 ± 0.22 
18.33 ± 0.20 
18.14 ± 0.19 
18.37 ± 0.20 
19.47 ± 0.22 
19.44 ± 0.23 
18.59 ± 0.22 
18.51 ± 0.22 
18.53 ± 0.22 
19.87 ± 0.28 
19.77 ± 0.27 
18.59 ± 0.23 
18.74 ± 0.25 
18.59 ± 0.23 
19.00 ± 0.29 
19.73 ± 0.25 
18.71 ± 0.24 
18.97 ± 0.28 
18.99 ± 0.28 

19.16 ± 0.30 
19.02 ± 0.28 



0.144 



+0.028 
-0.023 

0.082+^11 
0.076™ 
070+ 016 

n 100+0.025 
u - lzz -0.021 

u ' wo -0 .023 
3+0.024 
*-0.020 
,+0.014 
-0.011 
3+0.015 
-0.012 

0Qfi+ 022 

u.uyo_ 018 
o in4+ ' 023 

+0.023 
-0.019 
+0.013 
-0.096 
+0.013 
0.010 
+0.023 
0.018 



0.1193 
0.061^ 
0.063 4 



-0.018 
+0.020 



0.102 
0.042 
0.047; 
0.097 

0.084i°;°f 7 
0.096 
0.066 
0.048S 
0.086i8;8™ 

0.0681HS 

0-067±g;g?g 

n57+ 018 
u - uo ' -0.014 

065+ 019 
U.UDO_ 015 



/ 

V 
V 
V 

I 
I 
I 

V 
V 

I 
I 
I 

V 
V 

I 
I 
I 
I 
V 
I 
I 
I 
I 
I 



f Extinction correction: AY = 0.635, A{ = 0.371 X Flux density 



Table 7: MOA data table used in the V and I bands. Time is given since the GBM trigger. The 
exposure for each frame is 30 s. 



Days Filter Mag (Vega) Mag (AB) Flux i>Jy)f Ext. (Mag)} Exp. time (s) 



2.7404 cf 25.31 ±0.24 25.25 ± 0.24 0.54±g;if 0.68 4500 

2.7404 r' 23.66 ±0.10 23.84 ±0.10 1.63j*g;^ 0.47 4500 

2.7404 i! 23.19 ±0.11 23.60 ±0.11 1-82+°;^ 0.35 4500 

2.7404 z' 22.59 ±0.10 23.13 ±0.10 2.58+°;|> 0.26 4500 

2.7404 J 21.66 ± 0.29 22.59 ±0.29 3.79±J;|| 0.15 3600 

2.7404 H 20.92 ±0.28 22.31 ± 0.28 4.671};^ 0.10 3600 

2.7404 K > 18.2 > 20.1 < 35.5 0.06 3600 



f Flux density, corrected for Galactic extinction. 

t Th e Galactic extinction correction along the line of sight for E(B — V) 
19981 ) using a CCM Milky- Way extinction law ((Cardelli et al.lll989l ). 



0.175 (ISchlegel et al 



Table 8: GROND Data table; seeing during the observation ranged from l'/2 to l'/5 depending on 
the band. 



